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FOREWORD 


The  work  reported  herein  concerns  the  cushion  pressure 
characteristics  for  the  Air  Cushion  Landing  System  under  the 
conditions  of  stiff-operation,  and  was  performed  under  United  States 
Air  Force  Contract  No.  AF  33(615 )-69-C-1001  (Project  1369)*  with 
The  Ohio  State  University  Research  Foundation,  Columbus,  Ohio. 

The  content  of  this  part  of  the  report  deals  with  the 
theoretical  calculations.  The  second  part  of  the  report  will  be 
concerned  with  the  experimental  verification  of  the  theory  developed 
herein.  This  phase  of  the  work  was  carried  out  from  July  1969 
through  June  1970  and  the  report  was  released  by  the  author  November 
1970. 


This  report  has  been  reviewed  and  is  approved 


KENHERL1  H.  DIOGES 


Chief,  Mechanical  Branch 
Vehicle  Equipment  Division 
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ABSTRACT 


This  report  studies  the  performance  characteristics 
of  the  air-cushioned  landing  and  take-off  system  for  aircraft 
during  its  stiff-operation  mode.  Stiff-operation  is  obtained 
during  the  early  stage  of  the  take-off  period  or  the  later 
stage  of  landing.  Its  chief  feature  is  that  the  pneumatic 
supply  chamber  of  bleed  air  is  in  almost  parallel  configura¬ 
tion  with  the  ground.  The  supply  air  flows  vertically  down 
through  the  bleed  holes  and  is  then  deflected  outward. 

This  part  of  the  report  contains  the  theoretical 
treatment  of  the  problem.  The  results  are  in  the  form  of 
a  cushion  pressure  ratio  in  terms  of  the  supply  (trunk) 
pressure.  Analysis  was  performed  based  on  Incompressible 
viscous  theory.  The  second  part  of  the  report  to  be 
published  shortly  will  detail  the  experimental  results  in 
comparison  with  theory. 
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I  INTRODUCTION 


The  air-cushioned  landing  system  for  aircraft  is 
based  on  a  simple  fluid  dynamic  phenomenon  -  the  momentum 
principle.  It  has  the  all  important  advantage  of  being 
able  to  operate  over  terrains  other  than  paved  runways  or 
over  water.  Furthermore,  it  obviates  the  complicated 
retractable  gear  system. 

The  construction  consists  of  a  retractable 
pneumatic  bag  shaped  in  a  toroidal  manner  and  fastened  to 
the  fuselage.  The  pneumatic  ag  made  from  rubber-like 
elastic  material  is  provided  with  air  bleed-holes  and  can 
be  inflated  by  a  low-pressure  high-volume  air  source. 

During  operational  periods,  i.e.,  takeoff  or  landing,  the 
inflated  bag  forms  a  seal  between  the  supporting  base  area 
of  the  aircraft  and  the  ground.  In  this  manner,  a  cushion 
space  is  formed.  The  pressurized  air  is  forced  through  the 
distribution  holes  and  Is  deflected,  because  of  the  symmetry 
in  geometry  to  form  a  peripheral  ground  jet.  The  outward 
momentum  of  the  ground  jet  creates  a  positive  air  pressure 
In  the  cushion  space  or  the  base  region.  It  is  this 
positive  pressure  that  supports  the  aircraft  weight  during 
takeoff  or  landing  operation.  Figure  1  depicts  the 
schematics  of  such  a  system  together  with  the  current- 
nomenclature. 
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FIGURE  I  SCHEMATICS  OF  AIR-CUSHIONED 

LANDING  SYSTEM 


y 

I  DODrtl  1C  m  ATC 


FIGURE  2.  SIMPLIFIED  FLUID  DYNAMIC  MODEL 

FOR  STIFF  OPERATION 
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While  the  practicality  of  such  a  landing  system 


has  been  amply  demonstrated  in  successful  flight  tests 
by  the  Bell  Aerosystems  program,  further  research,  testing 
and  gathering  of  engineering  data  are  required  in  order 
to  obtain  more  precise  information  for  the  optimum  design 
of  such  a  system  for  various  modes  ^f  operation.  Some 
of  the  obvious  modes  of  operations  may  be  mentioned:  viz. 
during  landing  operation  the  attitude  of  the  aircraft  may 
necessitate  a  landing  with  a  fixed  pitch  angle  or  a  fixed 
roll  angle.  The  cushion  seal  is  no  longer  of  a  uniform 
height  around  the  periphery.  This  mode  of  operation  may 
considerably  alter  the  fluid  dynamic  behavior  near  the  base 
region  and  hence  the  maneuverability  of  the  aircraft. 
Alternatively  the  ground  Jet  distribution  may  be  influenced 
by  a  crosswind  which  again  may  change  the  handling  character¬ 
istics  of  the  aircraft  in  question.  Aside  from  these 
factors  it  is  to  be  noted  that  the  locations,  spacings  and 
sizes  of  the  bleed-holes  have  not  yet  been  optimized 
although  some  preliminary  work  was  reported  by  Earl  in  the 
report  by  Lell  [1]#.  Furthermore,  it  was  observed  that  the 
pneumatic  bag  experienced  some  breathing-mode  vibrations 
during  its  inflation  stage  which  is  indicative  of  the 
complexities  of  the  mechanics  of  deformation  of  such  a 


•Numbers  in  brackets  are  reference  numbers. 
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toroidal  elastic  enclosure.  The  exact  shape  of  the  cushion 
bag  under  pressure  may  be  of  importance  in  that  the  cushion 
pressure  distribution  in  the  base  area  may  be  greatly 
influenced  by  the  deformed  shape.  Similarly,  a  number  of 
design  considerations  dealing  with  the  controlled  bleed- 
flow  distribution  in  various  portions  of  the  pneumatic  bag 
may  be  mentioned.  The  controlled  distribution  may  be 
accomplished  by  compartmentalizing  the  pneumatic  bag;  each 
compartment  having  its  own  pressure.  This  scheme  can  be 
accomplished  by  a  proper  valving  and  ducting  and  can  be 
used  to  control  and  direct  the  aircraft  attitude.  The 
preceding-mentioned  phenomena  are  noted  herein  to  indicate 
the  need  for  further  studies  and  refinements  and  will  be 
more  fully  elaborated  upon  in  a  later  report. 

The  present  study  described  in  this  report  is, 
however,  concerned  with  a  particular  mode  of  operation 
known  as  the  stiff-operation  which  takes  place  when  the 
pneumatic  bag  is  very  close  to  the  ground.  The  air  seal 
formed  becomes  a  long  (in  the  air  flow  direction)  passage 
with  very  small  height  (clearance  to  the  ground).  Figure  1 
shows  the  schematics  of  the  various  major  components  comprising 
the  air-cushion  system.  For  the  stiff-operation  mode,  the 
bleed  air  is  directed  along  the  small  passage  outward 
leaving  the  cushion  space  relatively  undisturbed.  Consequently, 
the  fluid  dynamic  behavior  can  be  very  well  simulated  by 

k 


assuming  the  cushion  space  to  be  a  stationary  area  and  the 
two  air  passages  can  be  brought  together  as  shown  in  Figure 
2  in  which  the  ground  is  replaced  by  a  solid  plane  on  the 
x-axis  and  the  source  of  air  injection  is  replaced  by  a 
porous  plate  at  a  fixed  distance  H  above  the  ground.  The 
plate  has  a  length  of  2L.  The  injected  air  is  divided 
along  the  y-axis.  The  pressure  at  either  end  of  the 
porous  plate  is  taken  to  be  atmospheric.  The  objective  of 
the  present  study  is  to  analyze  the  pressure  along  the 
passage.  Since  at  x  «  0  along  the  y-axis  the  flow  is  zero 
due  to  symmetry,  hence  the  pressure  there  shall  be  taken 
as  the  pressure  in  the  cushion  space  which  is  assumed  to 
have  no  motion.  The  injection  of  air  is  to  take  place  with 
a  uniform  rate  independent  of  x.  The  latter  assumption 
will  be  removed  later  in  this  report.  A  subsequent  report 
will  discuss  and  ana'yze  the  fluid  dynamics  when  the  former 
assumption  is  removed,  i.e.  there  is  induced  motion  of  the 
air  in  the  cushion  space. 

II.  ANALYSIS 

Employing  the  assumption  of  an  incompressible 
fluid  model  and  two-dimensional  flow  tne  equations  of  motion 
and  of  continuity  are* 

2 

uu  +  vu  ■  -p  /p  +  vV  u 

A  J  A 

2 

UVX  +  Wy  «  -Py/P  ♦  W  V 


u  +  v,  .■  0 
x  y 
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(1) 

(2) 

(3) 


The  boundary  conditions  associated  with  the 


preceding  equations  are: 


u  -  v  ■  0 

y  -  0 

(4a) 

u  ■  0 

y  -  H 

(4b) 

V  -  -Vw 

y  -  H 

(4c) 

III.  UNIFORM  INJECTION  CASE 

Since  the  flow  in  the  passage  is  caused  by  the 
injection  at  the  upper  boundary  the  flow  velocity  in  the 
x-directi :n  must  be  linearly  proportional  to  x,  the  distance 
from  the  line  of  symmetry,  for  a  constant  rate  of  injection, 
i.e.  v  ■  -v  at  y  ■  H.  Consequently,  the  stream  function 
f  can  be  written  as: 

f  -  vw  x  F(y )  (5) 

where  y  ■  y/K 

v  «  -vw  F  (6) 

u  -  (vwx/H)  F '  (7) 

Substitution  of  the  above  relations  into  equations  (1) 
and  (2)  yields: 

F'2  _  pp«  .  -(H2p  /pv2x)  -  F"  '/Re  (8) 

x  w 

FF’  -  -(H  pv/pv2)  -  FH/Re  (9) 

where  Re  •  (H  vw/v) 

Equation  (9)  states  that  p^  is  independent  of  x; 
hence  px  is  independent  of  y.  In  equation  (8)  the  first 
term  on  the  right-hand  side  is,  therefore.  Independent  of 
y  and  can  be  regarded  as  a  constant,  or 
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(9a) 


-(H2  px/pv2x)  *  C 

the  principal  equation  concerned  in  this  report  is  then 

p*2  ^pn  «  c  +  F*  "/Re  (10) 

the  boundary  conditions  are 

y(o)  -  p»(o)  -  p»(i)  =  1  -f(i)  »  o  (id 

In  equation  (10)  the  Reynolds  number  is  a  parameter 
to  be  arbitrarily  specified  and  there  are  four  boundary 
conditions  in  (11)  for  the  third-order  equation.  The  constant 
C  analogous  to  an  eigen-value  in  a  linear  differential 
equation  is  to  be  determined  as  well.  Its  value  is,  of 
course,  dependent  on  the  Reynolds  member. 

Low  Reynolds  Number  Case 

For  very  low  Reynolds  numbers  the  Inertia  terms 
on  the  left-hand  side  of  equation  (1)  can  be  ignored 
since  they  are  of  the  second-order  nature  and  the  equation 
reduces  to 

P" '  *  -C  Re  (12) 

with  a  solution 

F  -  -(2y3  -  3y2)  (C  Re/12)  (13) 

which  satisfied  the  first  three  boundary  conditions  of 
(11)  and  the  last  condition  is  fulfilled  by 

C  -  12/Re  (1H) 

The  velocity  profile  from  7  ■  0  to  y  «  1  is 

given  by 

P»  -  y  (1-y )  (C  Re/2)  (14a) 

which  is  symmetrical  about  the  mid-channel  position  y  •  1/2 


t 


thus  exhibiting  no  convective  effect  from  the  injection 
•side  y  =  1. 

The  velocity  profiles  for  very  low  Reynolds  number 
are  graphically  shown  in  figures  3,  4,  ...  for  increasing 
values  of  the  Reynolds  number.  The  profiles  exhibit 
symmetry  about  the  midpoint  of  the  flow  channel,  char¬ 
acterizing  the  absence  of  the  inertia  effect.  The  pressure 
distribution  along  the  channel  is  solely  caused  by  the 
frictional  effect  and  can  be  calculated  directly  from 
equations  (9a)  and  (14).  Equating  these  two  equations, 

-H2p„/(p  v2  x)  -  C  -  12/ (vH/ v) 

A  w  W 

Denoting  the  atmospheric  pressure  by  p  ,  the  pressure 

a 

distribution  is  obtaineu  as 

(P-Pa)  *  (6y  vw/H3)  (L2  -  x2) 

The  maximum  pressure  occurs  at  the  mid-channel  x  »  0 

where  the  velocity  is  zero.  In  an  approximate  sense  it  is 

also  the  pressure  in  the  cushion  space.  Denoting  U  by  p  , 

c 

the  value  is 

Pc  -  Pa  -  vw  L2/h3 

This  equation  may  be  used  either  to  calculate  the  cushion 
pressure  or  the  injection  velocity  v  . 

Arbitrary  Reynolds  Number  Case 

Studies  similar  to  the  present  problem  have  been 
undertaken  by  a  number  of  investigators  although  for  different 
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applications  from  that  contemplated  in  the  present  case. 

Yuan  and  Pinkelstein  [2]  studied  the  problem  of  flow  in  a 
circular  tube  with  porous  wall.  The  resulting  equation 
governing  the  velocity  distribution  is  quite  similar  to 

(10)  and  (11).  Their  method  consists  of  solving  for  the 
stream  function  in  an  ascending  series  of  Re  for  small 
Reynolds  number  or  series  of  (1/Re)  for  large  Reynolds 
numbers.  Because  of  the  series  method  of  solution  a  large 
number  of  terms  are  required  for  accurate  evaluation  of 
the  velocity  profiles.  Consequently,  Morduchow  [3]  devised 
an  integral  method  modeled  after  the  K<  m-Pohlhausen 
method  to  solve  the  problem  posed  in  [2].  The  method  is 
well  known  to  be  reliable  only  when  it  can  be  shown  to 
agree  with  some  known  exact  solutions.  A  more  rigorous 
treatment  of  the  similar  problem  has  recently  been  given 

by  White  et  al  [4].  The  last-mentioned  work  solved  the 
problem  by  an  infinite  series  in  terms  of  the  physical 
coordinates.  The  method  is  conceptually  exact  but  is  tedious 
to  execute. 

In  the  present  study  the  method  of  solution  of 
equations  (10)  and  (31)  are  made  exact  both  conceptually 
and  numerically.  The  execution  of  the  solution  turns  out 
to  be  quite  simple.  In  the  first  place  equations  (10)  and 

(11)  constitute  a  two  fixed-point  boundary  value  problem 
with  an  associated  eigen-value.  The  "eigen-value"  C  in 
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equation  (10)  together  with  the  arbitrary  parameter  Re 
appears  to  be  complicated  but  in  reality  they  simplify 
the  problem  as  can  be  seen  as  follows: 


Let  F{y)»e  £(yy) 

n  *  yy 

Equation  (10)  becomes 

32  y2(f'2  -  ffw)  «  C  +  8y3  f"*/ Re 
Putting  g2  y2  *  C  ■  Sy^/Re  the  following  results: 

(15a) 

8  =  C1//4/Re1/2 

(16a) 

Y  «  C1/4  Re1/2 

and  the  equation  for  f  now  reads 

(16b) 

f’2  -  ff"  ■  1  +  f" » 

(17) 

which  turns  out  to  be  the  well-known  Falkner-Skan 

However,  the  boundary  conditions  are: 

f(0)  *  f»(0)  «  0 
f(C1/4  Re1/2)  -  Re1/2  C1/l 4 

equation. 

f»(C1/4  Re1/2)  »  0 

(18) 

Now  the  behavior  ox’  the  function  f  subject  to  the  first 
two  boundary  conditions  of  (18)  is  well  known,  ^or  f "(0) 
>  1.2326,  f’  increase  from  zero  monotonically .  For 
f" (0)  <  1.23?6,  the  graph  of  f  *  n  exhibits  a  maximum 
with  f’  <  1  and  the  f’ -curve  crosses  the  n-value.  Taking 
advantage  of  these  known  properties,  equation  (17)  with 
(18)  as  its  boundary  conditions  may  be  solved  by  assuming 

b 

an  arbitrary  value  of  fw(0)  between  0  and  1.2326. 
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PERIPHERAL  JET  VELOCITY  ( u/u0) 

FIGURE  3  VELOCITY  DISTRIBUTION  ACROSS  CHANNEL- SLAL  (Re  =  0) 


INJECTION  SIDE 


oocop-ioio^ion  —  o 
(H/A)  ONDOU9  WOUJ  J.H9I3H 
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PERIPHERAL  JET  VELOCITY  (u/u0) 

FIGURE  4  VELOCITY  DISTRIBUTION  ACROSS  CHANNEL-SEAL  (Re«U6) 


INJECTION  SJDE 


PERIPHERAL  JET  VELOCITY  (u/uG) 

FIGURE  5  VELOCITY  DISTRIBUTION  ACROSS  CHANNEL- SEAL  CRe=  2.17) 


INJECTION  SIDE 


PERIPHERAL  JET  VELOCITY  (u/ua) 

FIGURE  7  VELOCITY  DISTRIBUTION  ACROSS  CHANNEL- SEAL  (Res7.l6) 


INJECTION  SIDE 


PERIPHERAL  JET  VELOCITY  (u/ua) 

FIGURE  8  VELOCITY  DISTRIBUTION  ACROSS  CHANNEL- SEAL  (Re 


FIGURE  9  VELOCITY  DISTRIBUTION  ACROSS  CHANNEL-  SEAL  (Re  s  22.28) 


PERIPHERAL  JET  VELOCITY  <u/u0) 

FIGURE  IO  VELOCITY  DISTRIBUTION  ACROSS  CHANNEL-  SEAL  {  R«  «  32.62) 


INJECTION  SIDE 


PERIPHERAL  JET  VELOCITY  (u/u„) 

FIGURE  II  VELOCITY  DISTRIBUTION  ACROSS  CHANNEL- SEAL  (R®“  44.80) 


INJECTION  SIDE 


INJECTION  SIDE 


-  INJECTION  SIDE 


(H/A)ONnOUO  WOUi  1H9I3H 
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PERIPHERAL  JET  VELOCITY  (u/u0) 

FIGURE  18  VELOCITY  DISTRIBUTION  ACROSS  CHANNEL -SEAL  (R«  =  2602.03) 
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FIGURE  19  PRESSURE  GRADIENT  PARAMETER  vs.  REYNOLDS  NUMBER 
IN  SEAL™ CHANNEL .  C  =  ™(H2Px/^v£  X  ) 


Table  I 


Summary  of  Numerical  Results,  Uniform  Injection  Case 


f"(0) 

C 

Re 

F"(0) 

F”  (1) 

(C/C») 

.100 

90000 

•13333xlO-3 

6.0000 

-6.000 

1.0000 

.200 

5624.0 

.0021345 

6.0000 

-5.999 

1.0000 

.300 

1110.1 

.01083 

6.0050 

-5.994 

1.0021 

.1100 

350.59 

.034456 

6.0157 

-5.981 

1.0066 

.500 

143.02 

.085284 

6.0389 

-5.954 

1.0165 

.600 

68.46 

.18143 

6.0829 

-5.903 

1.0351 

.700 

36.497 

.35124 

6.1602 

-5.616 

1.0683 

.800 

20.975 

.64428 

6.2934 

-5.672 

1.1261 

.850 

16.235 

.86480 

6.3932 

-5.569 

1.1700 

.900 

12.702 

1.1618 

6.5271 

-5.438 

1.2298 

.950 

10.023 

1.5727 

6.7111 

-5.270 

1.3136 

1.00Q 

7.9601 

7.1656 

6.9738 

-5.051 

1.4365 

1.050 

6.3491 

3.0815 

7.3724 

-4.762 

1.6304 

1.100 

5.0723 

4.6732 

8.0407 

-4.375 

—  ' 

1.120 

4.6405 

5.6978 

8.4533 

-4.184 

— 

1.11*0 

4.2430 

7.1594 

9.0176 

-3.970 

— 

1.160 

3.8777 

9.4212 

9 . 8370 

-3.732 

— 

1.180 

3.5363 

13-403 

11.140 

-3.472 

— 

1.190 

3.3725 

16.800 

12.138 

-3.334 

— 

1.200 

3.2100 

22.281 

13*584 

-3.190 

— 

1.205 

3.1280 

26.514 

14.594 

-3.114 

— 

1.210 

3.0444 

32.621 

51.928 

-3.036 

0.9886 

1.215 

2.9577 

42.196 

17.800 

-2.953 

0.9898 

1.216 

2.9400 

44.803 

18.274 

-2.936 

0.9902 

1.217 

2.9217 

47.743 

18.792 

-2.918 

0.9904 

1.218 

2.9033 

51.086 

19.362 

-2.900 

0.9908 

1.219 

2.8846 

54.920 

19.995 

-2.882 

0.9910 

1.220 

2.8655 

59.361 

20.702 

-2.863 

0.9916 

1.222 

2.8260 

70.754 

22.404 

-2.824 

0.9925 

1.221* 

2.7843 

87- -36 

24.670 

-2.783 

0.9936 

1.226 

2.7393 

114.21 

27.900 

-2.739 

0.9948 

1.228 

2.6890 

164.25 

33.048 

-2.689 

0.9961 

1.230 

2.6293 

291.28 

43.345 

-2.629 

0.9976 

1.231 

2.5920 

474.31 

54.767 

-2.592 

0.9986 

1.232 

2.5415 

1277.3 

88.628 

-2.542 

0.9993 

1.1322 

2.5272 

1933.1 

108.59 

-2.527 

0.9995 

1.2323 

2.5188 

2602.0 

125.68 

-2.519 

0.9997 

1.2321* 

2.5087 

3980.8 

155. 00 

-2.509 

0.9998 

’  232587  »2/4 


v*  is  the  asymptotic  value  of  C  as  Re+  0  or  Re 
respectively  by  equations  (14)  ancl  (40). 


given 


With  f(0)  ■  f*(0)  •  0,  straightforward  integration  of  the 
equation  will  bring  f*  to  zero  at,  say,  n  ■  n0  at  which  the 
value  of  f  is,  say,  D.  In  view  of  the  last  two  conditions 
of  (18/  it  is  obvious  that 

C1/1*  Re1/2  -  n 


Re1/2/C1/4  ■  D 


and,  of  course. 


C  -  (n  /DT 

0  (19a, b) 

Re  ■  n  D 
0 

The  only  inconvenience  in  the  above  scheme  is  that  the 
Reynolds  number  is  the  outp  t  instead  of  being  an  input 
and  may,  therefore,  be  a  little  awkward.  It  goes  without 
saying  that  f"(0)  ■  1.2326  corresponds  to  Re  ■  •.  The  value 
of  C  is  shown  later  to  approach  w2/4.  In  short,  for  each 
value  of  fM (0)  assumed  there  results  the  following:  Reynolds 
number,  the  pressure  gradient  parameter  C,  the  velocity 
profile  F'  and  the  transverse  pressure  profile  as  supplied 
by  equation  (9).  The  results  of  these  calculations  are 
shown  in  Table  1  where  the  values  of  F"(0)  and  F"(l)  are 
proportional  to  the  shear  stresses  on  the  lower  wall  and 
the  upper  wall  where  the  injection  takes  place.  The  re¬ 
duction  of  the  shear  on  the  injection  side  is,  of  course, 
expected.  Of  interest  are,  however,  the  asymptotic  value 
of  the  pressure  parameter  C  and  the  shear  on  the  injection 
side.  Of  equal  significance  is  the  fact  that  at  Reynolds 
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numbers  larger  than  20  there  is  a  rough  relationship 
expressed  by 

F"(0)~  Re1/2 

This  relationship  and  th  gradually  steep,  velocity  gradient 
for  increasing  Reynolds  number  cn  the  solid  wall  as  shown 
in  Figures  15  to  19  suggest  the  existence  of  a  thin  boundary 
layer  near  the  solid  wall,  i.e.  the  ground. 

The  pressure  gradient  parameter  C  defined  in 
equation  (9a)  is  plotted  in  Figure  20.  It  is  clear  that 
at  the  lower  end  of  the  Reynolds  number  scale  the  asymptotic 
formula  C  -  12/Re  holds  whereas  at  the  high  end  the  C-value 
tends  to  be  a  constant  which  is  shown  in  the  next  section 
to  be  w  /k.  The  utilization  of  the  C-values  lies  in  their 
relation  to  the  pressure  distribution  along  the  flow 
channel  to  be  discussed  later. 

Asymptotically  Large  Reynolds  Numbers  Case 

In  common  with  the  concept  of  the  boundary  layer 
theory  the  term  F"'/Re  in  equation  (10)  is  dropped  and  the 
resulting  equation  will  be  treated  as  the  first-order 
external  (to  the  boundary  layer)  solution.  The  resulting 
equation 

pt2  _  pp*  .  c  (2o) 

has  tv.o  solutions,  i.e. 

f  -  ±vr? 
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and 


F  *  A  cos  y  Vc/(A°  +  L°)  +  B  sin  y  Vc/(A2  +  b2) 

The  frist  solution  is,  of  course,  a  member  of 
the  second  by  letting  A  *  0  and  B  -+•  »•  This  choice  of 
A  and  B  leading  to  a  flat  profile  for  F*  is  not  consistent 
with  the  results  of  the  numerical  solution  and  is  hence 
rejected.  Now  in  the  second  solution  there  are  three 
arbitrary  constants  which  are  to  be  determined  by  satisfying 
only  three  conditions*  of  (11).  Since  on  the  surface 
y  *  0  there  is  a  boundary  layer  which  is  to  be  the  solution 
of  a  boundary  layer  equation  the  conditions  F'(0)  *  0  is,, 
therefore,  dropped,  the  resulting  solution  is 

F  -  sin  y\ZTT  (21) 

where Vc-*  */ 2.  The  value  of  C  appears  to  approach 
(ir2/4)  in  Table  I.  This  solution  (external)  gives  F*(0)  * 
tt/2  and  is,  therefore,  not  valid  near  y  *  0.  In  order 
to  construct  a  sequence  of  solutions  in  an  orderly  manner 
equation  (10)  is  solved  by  the  following  expansion  In  the 
external  region: 


C1/*  -  cj/l4  [l  +  Cx  e  +  C2  c2  ...j 

(22) 

C  -  »2/4  andc  *  1//We 

(23) 

0 

F  -  Fo(y)  +  cF1(y)  +  c2  F^y) 

(24) 

*The  problem  treated  by  White  et  al  [4]  is  the  one  in  which 
all  their  four  conditions  can  be  satisfied  by  three  arbitrary 
constants » 
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After  the  usual  procedure  there  results: 


F»2_pp,i*c 
o  oo  o 

2F«  F»  -  F  F"  -  F”  F. 
01  oi  ol 


4C  C 
1  o 


(25) 

(26) 


2  2 
Fl'  +  2Fo  F2  "  P?V'  "  F1F1  "  F'52  "  C0(4C2+6C^)  +Fq  * ( 27 ) 

In  the  boundary  layer  region,  l.e.  in  the  vicinity 

1/4  - 

of  Jf  ■  0  the  boundary  lav.  er  coordinate  ■  C  y/e  is 

u  o 

used  and  expand  the  function  F  as 


F  *  fn(n  '  +  €  f ( n0 ^  +  e^f«^n  )  +  •••] 

1/4  r  0  1  ,  '•  0  J  . 

fo  +  t(clf0tfl)  +  E‘  (f2+0lfl+C2f'o)  +  •••! 


eC 


F'(y)  -  C1/2  If*  +  e(C  f '  +  f 1 )  +  ...1 
oLo  lol  J 

The  composite  function  f  ,  f  ,  ...  are  determined  by: 


(28) 

(28a) 


f -  f 
o  o 


,2  _ 


f  f"  -  1 

0  0 


(29) 


f”'+f0f”  -  2rf»  +  f«  f1  •  ^•4c1  -  2C1f0f”+  2c1r;2  -  c1f”«  (30) 

It  is  to  be  ncoed  that  equation  (29)  although  the 
same  fora  as  equation  (17)  is  in  reality  different  in  that 
the  former  is  the  boundary-layer  equation  subjected  to  the 
standard  two-point  boundary-value  problem  whereas  the  latter 
is  valid  for  the  entire  channel  and  constitutes  an  eigen¬ 
value  problem.  Beginning  with  (25)  one  has 


F0  -  sin  y  (31) 

where  the  parameter  CQ  is  taken  to  be  n/2  in  order  to 
satisfy  the  condition  F(l)  ■  1,  hence 


33 


c  »  *2/4  (32) 

o 

Near  the  solid  wall  $  *■  0  equation  (31;  gives  a  finite 
velocity  F^(0)  86  (n/2).  Thi3  is  to  be  remedied  by  allowing 
a  boundary  layer  to  exist  near  y  *  0  of  th*  thickness^  in 
which  the  velocity  P'  is  to  vary  from  0  to  (n/2)  at  the 
edge  of  the  boundary  layer  6.  This  behavior  is  to  be 
provided  by  the  first  term  f  in  the  expansion  near  y  *  0‘ 
as  shown  in  equation  (28).  The  governing  equation  (29) 
for  f  is,  therefore,  subjected  tc  the  boundary  conditions 

l/ii“ 

fo(0)  -  0;  f£(0)  *  0;  and  f^(C^  6/e)  *  1 

At  asymptotically  large  .  nolds  number,  i.e.  e  -*0,  the 
last  condition  is  replaced  by  f”(“)  ■  1.  The  solution  to 
equation  (29)  with  these  boundary  conditions  is,  of  course, 
the  Falkner-Skan  function  with  f" (0)  *  1.2326  and  the  value 
of  f  at  large  arguments  is  asymptotically  given  by 

f  =  n  -  0.65 

o  o 

This  composite  nature  of  the  function  P  comprised  of 

equation  (31)  in  the  main  bulk  of  the  channel  and  the 

Falkner-Skan  velocity  distribution  near  the  solid  wall 

is  illustrated  in  Figure  21.  In  the  boundary  layer  region 

y  *d  the  dashed  curve  e  Cly^d  f  /dy  *  C1//2  f  replaces  F*. 

q  o  o  o  o 
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The  function  F  Itself  In  this  region  is  replaced  by  e  C^'  " 

1  /i\ 

However,  this  replacement  i.e.  e  C  '  f  does  not  match 
*  *  o  o 

the  externa]  stream  function  at  the  "junction"  y  =  5.  The 
disparity  i3  given  by 

T  /i4  1/4  _ 

AF  *  F^  (6)  -  e  C"'  f  (C  6/e) 

O  o  o  o 

Under  the  condition  of  6+  0  and  6/e-*«  the  disparity  is 
obtained  as 

AF  *  •''cT  6-  e  C1/li  «A  -  0.65!  -  0.65  e  C1/i4 

o  o  L  o  J  o 

Physically  this  quantity  represents  the  transverse  displace 

ment  velocity  due  to  the  boundary  layer  near  the  solid  wall 

The  mismatch  of  the  stream  function  F  at  y  =  6 

is  now  to  be  corrected  in  the  outer  region  of  eF^,  the 

second  term  of  equation  (2'D,  such  that  the  value  of  F 

1/4 

+  eF^  is  equal  to  that  of  e  C0  fQ  at  y  8  i.  This  is 
tantamount  to  the  depression  of  the  F-value,  to  be  accompli 
shed  by  cF^,  by  an  amount  equal  to  AF.  Consequently,  the 
conditions  of  F^  are: 

1/4 

F, (1)  =  0;  F ' ( 1)  =  0;  and  F,<6)  *  >0.65  C 
i  1  j-  o 

For  small  6  the  last  condition  is  replaced  by 

F1(0)  =  -0.65  cj/4  (3 

The  solution  of  equation  (26)  with  the  first  two  boundary 


conditions  is 


?1(y)  =  -2C,  f(l-y)  Vc;  cos  *<|  (3*0 

and  for  (33)  to  hold  one  has: 

C1  =  0.65/(20^)  (35) 

Now  the  process  repeats  itself  i.e.  F^(y)  gives  a  slip 
velocity  at  y  =  0  as 

e  F*  (0)  =  0.65  C1/if  e  (36) 

1  o 

which  to  be  connected  by  the  second  term  in  the  inner 
expansion  of  equation  (28),  i.e., 

e  C1/2  |c  f*  +  f'l  =0.65  C1/1<  e 
o  L  1  o  lj  c 

at  n  .  Hence,  the  boundary  condition  of  f  at  n-*-®  is 
0  1 
simply  f|(  )  =  C^.  The  other  two  conditions  are,  of 

course,  f’(0)  =  0  and  f  (0)  =  0  and  the  solution  f 

satisfying  equation  (30)  turns  out  to  be 

fl  *  °1  "o  r0  (37) 

The  outer  stream  function  and  the  inner  stream  function 
at  y  =  differ  by 

AF  =  { f  0  ( 6 )  +  e  F^(6)|  -  eC^  jf^®)  +e  jc^f^C®)  + 

=  0.652  e2/2 

which  is  used  to  determine  the  boundary  condition  for  F^ 
at  y  =  6,  i.e. 

F2(0)  =  -0.652/2 
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The  function  F2  from  equation  (27)  is  found  to  be 


,(y)  -  -(C1  TT2/2)  [(1-y)  sin  (y/2)]  +  (ir/2)  [4^  • 

2C2  -  (l-y)cos(*y/2)  -  (ir/2)cos(y/2)^Try^ 


(  tt/2  ) 


£cos  z(sin  z  -  z  cos  z)/sin3zj  dz  +  (?r/4) 

|Vin(iry/2)-(  iry/2)cos  ( ir/2 )j  j^(  cot  (ity/2 )/ sin(  uy/2 )  )+log 

i(ry/4)]  (38) 


tan  i 


The  condition  of  P2(0)  =  -0.65  /2  is  used  to  determine  the 


coefficient  C2  as 


2C 


n/2 

2  =  2C^  -  3C2  +  (0.652/tt)  +  \  ^cos  z  (sin  z  -  z  cos  z) 

0 

/sin^zj  dz 


or 


C2  =  0.4337  (39) 

Up  to  the  second  order  of  e  ,  the  pressure-gradient  constant 
C  may  be  expressed  as 

C  =  (tt2/4)(1  +  1.0373  e  +  2.1682  e2  ...)  (40) 

and  the  shear  stresses  at  y  «  0  and  y  =  1  are  given  by  the 
following  asymptotic  formulas 


F"(Q)  =  (C3/Ve  ) 

0 

[l  +  ?C^  £  +  .  .  .] 

(41) 

F"(1)  =  U2/4)  [: 

L  -  4C^  e  +  O(e^)  +  ...J 

(42) 
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The  pressure  differential  across  the  channel  is  given 
directly  as 

<Vi '  Vo5  =  ‘  pv«/2  (i,3) 

The  pressure  gradient  C,  the  shear  stresses  P"(0)  and  F” (1) 
and  other  pertinent  quantities  are  shown  in  Table  1. 

Figure  (20)  shows  the  variation  of  C  vs.  the 
injection  Reynolds  number.  From  an  engineering  calculation 
viewpoint  the  pressure  distribution  is  given  by 

P  -  C  [pv^/2]  (L2  -  x2)/H2  +  pa  (4i*) 

where  pa  is  the  ambient  pressure  and  the  use  of  this  relation 
implies  that  the  pressure  difference  across  the  channel 
width  as  given  by  equation  (43)  is  negligible.  At  x  *  0, 
the  peak  pressure  is  given  by  equation  (44)  with 

Pc  -  Pa  *  0  [pyjj/2]  (L2/H2)  (45) 

where  pc  now  denotes  the  cushion  pressure.  It  is  not 
possible  or  meaningful  in  the  present  case  to  construct 
the  pressure  parameter 

(pc  -  pa>/(pt  -  pa> 

where  pt  is  the  total  pressure  in  the  trunk.  The  reason 
is  that  the  injection  rate,  i.e.,  the  rate  of  air  flow 
through  the  trunk  distribution-holes  is  a  constant  is  only 
true  at  the  choked  condition.  However,  if  one  tentatively 
places  the  total  (trunk)  pressure  as 
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Pt  ■  P0  *  1/2  0  v» 


then  the  pressure  coefficient 

<pc-pa)/(pt-pa)  “  C(L2/H2)/  £l  +  C(L2/H2)]  ...  (46) 

which  is  seen  to  be  dependent  on  the  Reynolds  number. 

Under  the  almost  contact  condition  between  the  ground  and 
the  trunk  LVH  ,  and  the  pressure  coefficient  becomes 
one. 


IV.  VARIABLE-INJECTION  PLOW  CASE 

The  analysis  presented  in  Section  III  is  valid 
when  the  rate  of  the  bleed  air-flow  through  the  distribution 
holes  along  the  cushion  seal  is  uniform.  In  practice  this 
is  only  true  when  the  trunk  pressure  is  at  a  very  high 
pressure  in  the  order  of  50  to  100  psig  so  that  the  flow 
through  the  orifices  is  at  a  choked  condition.  Since  this 
extreme  condition  is  rarely  encountered,  it  is,  therefore, 
desirable  to  have  a  method  available  to  account  for  the 
variable  flow  case  so  that  more  realistic  cushion  pressure 
can  be  estimated. 

Prom  a  practical  consideration  the  distribution 
of  the  bleed  flow  from  the  pneumatic  trunk  is  affected  by 
(1)  the  flow  resistance  across  the  ’•orifices,”  i.e.  the 
orifice  coefficient  and  (2)  the  pressure  distribution  in 
the  flow  channel.  Thus  the  distribution  of  the  bleed  flow 
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along  the  channel  cannot  be  determined  until  the  pressure 
variation  in  the  channel  is  known;  the  latter  is  in  turn 
dependent  upon  the  flow  distribution  or  accumulation  from 
the  bleed-side  of  the  wall.  While  it  is  possible  to 
analyze  the  interaction  as  discussed  above  from  a  more 
exact  fluid  dynamic  viewpoint  it  is  more  expedient  to 
employ  the  concept  of  local  similarity  in  treating  the 
present  problem. 

The  local  similarity  rule  as  applied  to  the 

channel  flow  problem  is  interpreted  a c  follows:  referring 

to  Figure  2  the  local  pressure  gradient  p  is  determined 

x 

by  a  local  pressure  gradient  coefficient  C  which  is  in 
turn  determined  by  the  local  Reynolds  number  (vwH/v).  The 
velocity  distribution  (u/ua)  is  also  determined  by  the 
local  value  of  the  Reynolds  number.  Mathematically  this 
is  expressed  by 


-  rx 

F'(y)  l  (vw/H)dx 
*  0 

(47) 

-vw  F(y) 

(48) 

Equations  (47)  and  (48)  reduce  to  equations  (6)  and  (7) 

when  v,T  is  a  constant, 
w 

Substitution  of  equations  (47)  and  (48)  into 
equation  (1)  there  results 

F'2  -  FF"  -  -  [H2px/Pvwdx]  +  (v/Hvw)F"»  (49) 
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The  first  term  on  the  LHS  is  now  a  function  of  x.  This  is 
because  the  Reynolds  number  (v  H/v)  is  x-dependent  in  (49) 

W 

but  not  so  in  equation  (8)  for  the  uniform  injection  case. 

In  the  present  case  the  present  gradient  term  in  the  brackets 

of  equation  (49)  is  only  dependent  on  (Hv  /v).  The  fact 

w 

fhat  this  is  a  practically  acceptable  procedure  is  justified 
on  the  basis  that  the  average  value  of  the  convective  terms 
on  the  LHS  of  equation  (49)  is  nearly  a  constant.  The 
pressure  gradient  p  is  dependent  upon  the  friction  resistance 

A 

$ 

on  the  wall  denoted  by  F" '  in  equation  (49)  and  upon  the 
momentum  Increase.  The  momentum  increase  is  partly  dependent 
upon  the  accumulated  flow  and  partly  upon  the  shape  of  the 
velocity  distribution  in  the  channel.  Consequently,  the 
integrated  value  from  y  =*  0  to  y  *  1  of  the  convective 
terms  of  the  LHS  represents  the  influence  of  the  shape 
factor  or 

(F*2  -  FF")  dy  *  2  f1  F,2dy 

■'o 

From  Re  =  0  to  Re  ■*  this  shape  factor  changes  from  1.20 

o 

to  ( it  /8).  Since  this  change  is  minor  and  this  variation 
is  included  in  the  local  similarity  rule  the  procedure  is, 
therefore,  valid. 

Computation  Procedure 

With  the  pressure  gradient  parameter  C  defined  in 
Table  I  for  various  local  Reynolds  number,  the  pressure 
gradient  can  be  written  as 
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c 


d£ 

dx 


0VW  J  vM  dx 
0 


(5u; 


For  flow  with  appreciable  variations  in  density,  it  should 
be  inside  the  integral. 

It  is  assumed  that  the  following  data  are  known 
in  advance: 

L  «  length  of  channel j  ft 
H  «  ground  clearance,  ft 

p 

v  ■  air  kinematic  visco  ity,  ft  /sec 
Pt  *  trunk  pressure,  psfa 
p„  ■  ambient  pressure,  psfa  (2116.8  psfa  or 

u 

14.7  psia) 

Additionally,  appropriate  coefficients  of  discharge  for 
flow  through  the  distribution  holes  must  be  known.  The 
procedure  is  essentially  a  cut-and-trial  method.  It 
begins  with  an  assumed  value  of  the  pressure  at  x  *  o  which 
is  below  the  trunk  pressure  and  above  the  ambient  pressure. 

For  illustration  let  the  bleed  flow  be  determined 
from  the  simple  orifice  formula: 

Vw  "  [Cd  ArV'2g0(Pt-P)/c]  (51) 

where  A„  is  the  ratio  of  the  orifice  area  to  the  seal  area, 
r 

For  example  If  the  orifices  are  placed  3  inches  apart 
(center-to-eenter)  on  a  square-pitch  basis  and  the  orifice 
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diameter  is  2  inches,  then  Ar  =  £  22/4 J  /32  =  tt/9.  Equation 

(51)  gives  the  value  of  v  at  x  =  0  which  defines  the 

w 

Reynolds  number  (v  H/v)  at  x  *  0.  From  Table  I  the  pressure 
gradient  parameter  C  is  obtained.  Since  the  pressure 
gradient  itself  is  zero  at  x  *  0,  the  pressure  at  x  =  b 
may  be  calculated  as 

P  <x-»)  -  p  (x-0)  -  t»\C/H2)  t?n  (52) 

The  next  step  is  now  to  compute  vw  at  x  *  4  by  means  of 
equation  (51)  and  to  compute  the  pressure  gradient  from 
equation  (50A  after  Tab^e  I  is  used  to  obtain  C.  In  this 
way  the  pressure  at  progressively  large  x  can  be  numerically 
evaluated  until  at  x  *  L  the  claculated  pressure  must  match 
the  ambient  pressure.  If  the  calculated  pressure  is  higher 
than  the  ambient  pressure  then  the  assumed  pressure  at 
x  =*  0  is  too  large  and  a  lower  value  must  be  tried. 

A  large  number  of  computer  runs  were  performed. 

The  general  ranges  of  the  parameters  were  as  follows:  trunk 
pressure  varied  from  0.5  to  6.0  psig  at  1.0  psi  interval; 
the  C  (coefficient  of  discharge  times  the  area  ratio) 
varied  from  0.003  to  0.27  at  0.003  interval;  the  ground 
clearance  height  H  varied  from  0.0*1  to  0.16  ft  at  0.0*1  ft 
Interval  and  the  clearance-length  ratio  (H/L)  varied  from 
0.05  to  1.0  ut  0.05  interval. 
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Figures  21r-aJb?c1  contain  the  ty  .cal  pressure 
variation  along  the  channel-seal  for  an  (H/L)  ratio  of 
0.05,  0.10,  and  0.15  respectively.  The  abscissa  is  the 
dimensionless  distance  (x/L)  from  the  cushion-space  and  L 
is  the  length  of  the  channel-seal.  The  ordinate  shows  the 
non-dimensional  local  pressure  ratio  (P-P  )/(P^-Pa) .  The 
parameters  are  the  values  of  the  CdAr.  For  large  values  of 
C^A  the  pressure  distribution  is  within  the  flo*  channel 
almost  constant  indicating  the  phenomenon  that  air  injection 
takes  place  only  near  the  exit  (ambient)  side.  It  is 
significant  to  call  attention  to  tne  fact  that  these  non- 
dimensional  pressure  distribution  is  virtually  independent 
O'"  the  trunk  pressure  at  least  from  0.5  to  6.0  psig  and  the 
ground  clearance  H  from  0.0^  to  0.15  ft-  Over  these  ranges 
of  the  above-cited  parameters  the  curves  are  identical  and 
Indistinguishable  from  one  another.  Numerical  differences 
however  do  exist  but  only  in  the  third  significant  digit. 

Each  curve  in  Figures  21~a,b,c,  was  used  to  obtain 

the  average  pressure  within  the  channel  seal.  This  averaged 

pressure  when  multiplied  by  the  area  of  the  seal  gives 

additional  lift  to  the  vehicle  which  may  be  of  an  appreciable 

magnitude.  The  parameter  shown  on  the  ordinate  of  Figure 

22-a,b,  is  (P„  -P  )/(P  -P  )  with  the  cushion  gage  pressure 
av  a  c  a 

as  the  reference;  P_„  is  the  average  channel  pressure. 
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It  is  to  be  noted  that  Figures  22-a,  and  b  are  for  two  different 
trunk  pressures  of  6.0  and  2.0  psig  and  these  two  figures 
virtually  duplicates  each  other. 

From  the  pressure  distribution  like  those  in 
Figure  21-a,  the  value  at  (X/L)  =  0  is  the  cushion  pressure 
ratio  based  on  the  assump:ion  of  no  air  moisture  in  the 
cushion-air  space.  This  cushion  pressure  ratio  then  depends 
on  (H/L)  and  the  C^r-value.  For  example  at  H/L  =  .05  and 
*  *^3.  . 0 5 •  •  •  the  cushion  pressure  ratios  are  0.549, 
0.79C,....  The  latter  values  are  tabulated  in  Table  II  and 
plotted  in  Figures  23-a,b,  for  the  trunk  pressure  of  6.0 
psig  and  2.0  psig.  These  two  figures  are  almost  identical 
indicating  its  insensitivity  to  the  trunk  pressure.  The 
last  two  figures  are  particulary  Important  in  that  the  overall 
performance  of  the  cushion-lift  can  be  directly  obtained. 

As  an  illustration,  consider  the  example: 
lrunk  pressure  *  4.0  psig 
Cushion-spcce  area  *  200  sq.  ft. 

Channel-seal  length  *  1  f t . 

Vehicle  weight  ■  100,000  lb. 

C  d^r  -  *21 

The  last  value  can  be  estimated  from  the  geometrical  spacing 
of  the  faTeed  holes  located  In  the  cushion-bag.  The  coefficient 
of  discharge  is  usually  in  the  ~ange  of  0.6  corresponding  to 
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the  orifice  discharge  coefficient.  A  more  accurate  estimate 
will  te  reported  in  Part  II  of  thi3  report  to  be  published 
shortly . 

From  the  preceding  in'"  rmation  it  is  possible  to 
calculate  the  necessary  cushion  pressure  for  take-off 
operation. 

Pc  *  100,000/(200  x  144)  =  3 . 47  psig 
The  cushion  pressure  ratio  is  therefore  (3.47/4.0)  =  0.868. 

From  Figure  23a  or  23b,  the  curve  with  Qj Ar  =  0.21  gives  a 
value  of  (H/L)  «*  0.20.  For  L  -  1  ft,  the  ground  clearance 
height  H  =  0.20  ft  or  2.4  inches.  Of  course  the  information 
in  Figure  23  can  be  used  to  determine  the  proper  CdAr  if 
H  is  to  be  specified.  This  specifies  the  proper  spaeings 
of  the  bleed  holes-  and  subsequently  the  rate  of  discharge 
of  bleed  air. 

A  general  computer  program  is  attached  in  Appendix. 

The  computer  program  requires  the  following  input  values: 

CdAr;  trunk  pressure,  ambient  pressure,  both  in  psia; 

seal  length  in  ft.;  ground  clearance  H  in  ft;  and  the  kinematic 

viscosity  of  air  (% 0.00015  ft^/sec).  These  data  are  to  be 

read  by  the  computer  program  on  a  single  punch  data  card 

(format  6F10.2)  for  each  configuration.  After  one  run  is 

completed  the  program  will  read  the  next  card  which  may 

contain  different  parametric  inputs  for  the  next  run.  The 

program  continues  in  this  manner  until  it  reads  a  blank 

card  which  must  be  attached  next  to  the  last  data  card. 
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The  output  of  the  program  for  each  run  consists 
of  two  parts:  the  first  part  is  a  tabulation  of  the  local 
pressure  ratio  (P<-Pa)/(Pt*rPa) ,  whe  local  injection  velocity 
(Prorated  over  the  entire  seal  area)  in  ft  per  sec,  the 
local  pressure  gradient  constant  C,  and  the  local  Reynolds 
number  at  about  39  stations  of  increasing  (x  /L)  values  where 
x  is  the  distance  from  the  cushiony-space.  The  second  part 
gives  the  average  pressure  (gage)  in  the  seal  area  in  its 
ratio  to  the  gage  trunk  pressure,  the  actual  injection  flow 
rate  in  cfs  per  ft  of  depth  of  the  cushion  space. 
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FIGURE  2! a.  Pressure  Variation  Along  Channel 
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FIGURE  22  a.  Average  Seal  Pressure  As  A  Function 
Of  Area  Ratio  And  Ground- Clearance 
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FIGURE  22  b.  Average  Seal  Pressure  As  A  Function 
Of  Area  Ratio  And  Ground-Clearance 
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FIGURE  23a  Cushion  Pressure  Ratios  As  A  Function 
Of  Area^Ratio  And  Ground- Clearance 
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FIGURE  23  b.  Cushion  Pressure  Ratios  As  A  Function 
Of  Area-Rotio  And  Ground- Clearance 
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TABLE  II 

COMPUTER  CUSHION  PRESSURE  COEFFICIENTS 
BASED  ON  VISCOUS  THEORY 
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.0022  0.0049  0.0088  0.0136  0.0196  0.0267 


TABLE  II  (cont. 


03^5  0.0532 
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APPENDIX: 


FORTRAN  COMPUTER  PROGRAM 
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c****** instruct  ion 

C******COEFD  IS  PRODUCT  OF  ORIFICE  COEF  AND  A*EA  RATIO 
C******PTRNK  IS  THE  TRUNK  PRESSURE  IN  PSI  ABSOLUTE 
C******PAMB  IS  ATMOSPHERIC  PRESSURE  IN  PSI  ABSOLUTE 
C******EL  IS  THE  SEAL  CHANNEL  LENTH  IN  FEET 
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